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The many-body Green's function theory with the random-phase approximation is applied to the 
study of easy-plane spin-1/2 ferromagnets in an in-plane magnetic field. We demonstrate that the 



usual procedure, in which only the three Green's functions {(Sf ; 5* )} (a* : 



-,z) are used, yields 



unreasonable results in this case. Then the problem is discussed in more detail by considering all 
combinations of Green's functions. We can derive one more equation, which cannot be obtained 
by using only the set of the above three Green's functions, and point out that the two equations 
contradict each other if one demands that the identities of the spin operators are exactly satisfied. 
We discuss the cause of the contradiction and attempt to improve the method in a self consistent 
way. In our procedure, the effect of the anisotropy can be appropriately taken into account, and the 
results are in good agreement with the quantum Monte Carlo calculations. 

PACS numbers: 75.10.Jm, 75.70.Ak, 75.30.Gw, 75.40.Mg. 
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I. INTRODUCTION 

Many-body Green's function theory is a powerful tool 
for theoretical studies of spin systems. In this formalism, 
a certain decoupling approximation is generally required 
to terminate the infinite hierarchy of equations of mo- 
tion for higher-order Green's functions. The first (low- 
est) order decoupling scheme introduced by Tyablikov, 1 
which is called the random-phase approximation (RPA) 
or "Tyablikov decoupling," is a very simple yet effective 
way to perform this operation. Many previous authors 
have attempted to go beyond Tyablikov's method and 
generated a variety of decoupling procedures such as the 
Callen decoupling, 2 Tahir-Kheli's theory, 3 the modified 
versions of the Callen decoupling^ and Oguchi's varia- 
tional theory^ Nevertheless, it is known that the RPA de- 
coupling is still the simplest and most reliable first order 
approximation [see, for example, the comparison between 
the RPA, the Callen decoupling, and the quantum Monte 
Carlo (QMC) calculations in Ref. |7(. More quantitative 
results may be obtained by applying the second-order 
Green's function theory originally proposed by Kondo 
and Yamaji 8 - for the one-dimensional isotropic Heisen- 
berg model. In this theory, the hierarchy of equations 
of motion is terminated at the second step with intro- 
ducing vertex parameters. Shimahara and Takada- ap- 
plied this theory to the two-dimensional case, and Junger 
et alJ£ and Antsygina et alJ^- extended it to the case 
where a uniform external field is applied. Several other 
methods 7 - ' 12 : 13 have also been attempted. However, some 
problems, e.g., how to determine the vertex parameters, 
are still under discussion, and thus these formalisms have 
not been fully established yet. Therefore, the simple RPA 
decoupling scheme is still often used for analytical studies 
of complicated systems. 

Field-induced phenomena in spin systems have been 
the focus of both theoretical and experimental studies. 
For example, in uniaxially anisotropic Heisenberg anti- 



ferromagnets, the magnetic field applied along the easy- 
axis induces multicritical behavior at the triple point 
of the antiferromagnetic, spin-flop, and paramagnetic 
phaseSi 14 : 15 ' 16 : 17 ] 18 The spin reorientation transition in- 
duced by a transverse magnetic field has also attracted 
considerable interest for both ferromagnetic and an- 
tiferromagnetic casesi^^^^^^^^ 2 ^* 2 ^^ 
Some of them were studied by using Green's function 
formalism with the RPA decoupling. However, in this 
paper, we point out that one should pay careful attention 
in applying the RPA decoupling scheme to such a compli- 
cated system, which has an anisotropy in the plane per- 
pendicular to the direction of the magnetization. As an 
example of such systems, we consider spin-1/2 ferromag- 
nets with an easy-plane exchange anisotropy. This sys- 
tem can be described by a comparatively simple model, 
and thus we can clarify the issues. 

This paper is organized as follows. First, in Sec. [Til 
we introduce the model Hamiltonian considered in this 
paper. In Sec. IIIH we demonstrate the application of the 
usual procedure, in which only a restricted set of Green's 
functions is used. In Sec. IIV1 we show that one more 
equation, which cannot be obtained by the above pro- 
cedure, can be derived by considering all combinations 
of Green's functions. Moreover, we point out that the 
two equations contradict each other if one demands that 
the identities of the spin operators are exactly satisfied. 
In Sec. [V] the cause of the problems found in the pre- 
vious section are discussed, and we attempt to improve 
the method in a self consistent way. The obtained results 
are compared with the QMC data. Finally, a summary 
is presented in Sec. IVII 



II. MODEL HAMILTONIAN 

The many-body Green's function theory for spin sys- 
tems, which is briefly reviewed in the Appendix, has been 
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developed by many authors We consider the 
application of this theory to the study of easy-plane ferro- 
magnets in an in-plane magnetic field. The Hamiltonian 
of this system is given by 

H = \ E MsfSf+sisv+bstsfi-hY, sf, (i) 



and 



three Green's functions Gfj = {{S^;Sj)) v , G^ — 
((Sf\Sj)) v , and G^ = ((S-:S~)) rr The equations 
of motion for these Green's functions are given by 



W G£- M = {{S?,Sr] n ) + m,B\;Sr)) n>u , (4) 



with 



I s ? >H] = Tj2 Jik ( S * S t ~ AS t S k) T hSf, (5) 



Sj ] — i£fj,u\Sj Sij , 



(2) 



and 



where Si = (Sf,Sf,Sf) is the usual spin operator at 
site i and J+j is the exchange interaction strength be- 
tween site i and site j. We take into account only the 
nearest- neighbour coupling, i.e., Jy = J > if i and j 
are nearest-neighbour sites and Jy = otherwise. In 
the spin-1/2 case, the single-ion anisotropy energy, e.g., 
— K2,i{Sf ) 2 , is constant and does not affect our re- 
sults. The behavior of the system depends strongly on 
the value of the anisotropy parameter A: for < A < 1, 
the spins prefer to lie in the xy plane (easy plane) while 
for A > 1, the spins tend to align along the z axis (easy 
axis). In this paper, we focus on the easy-plane case. 
The magnetic field h in Eq. ([1]) is applied perpendicular 
to the hard (z) axis. 

Let us calculate the magnetic properties of the model 
given by Eq. (JT|) . Since the model is expressed in terms of 
spin operators Si, we have many choices of the operators 
A and B in Eq. (|A1[) . First, the calculations with the 
usual set of Green's functions are shown in the next sec- 
tion. Then, considering all choices of Green's functions, 
we discuss the problem in more detail in the subsequent 
section. 



[Sf,H] = ~J2 J *( S t S k SrS+)~(S+-Sr). (6) 



The right-hand sides of the above equations include 
higher-order Green's functions, e.g., ({SfS^',S7}) v ,ui- In 
order to close the system of equations, we adopt the 
generalized version of the so-called RPA or Tyablikov 
decoupling, 1 

((S£ S%; Sj)) VjU1 rj (Si)((S^;Sj)) riiU +(S^)((Si;Sj)) V}U i 
-(m)(SZ);S})) n , u . (7) 
Performing the Fourier transformations given by 

C uv _ (JH» -iq-(Ri-Rj) (o\ 



and 



(9) 



where N is the number of lattice sites, we now rewrite 
Eq. ([4]) in a matrix form: 



III. USUAL CHOICE OF GREEN'S FUNCTIONS 

In the case of easy-plane anisotropy (0 < A < 1), only 
the component of the magnetization parallel to the in- 
plane magnetic field has nonzero value, i.e, (Sf ) = m 
and (Sf) = (Sf) = 0. In order to study the magnetic 
properties of spin-S* systems, the set of Green's functions 
, G;. ' m ™, and Gf~' m ™ has been often used in 



G+ 



many previous studies. Here we denote 



ij,V 



= ((s>t;(s*r(s-r)) 



(3) 



where m > and n > 1 are positive integers. S^ = 
Sf ± iSf are the usual spin raising and lowering opera- 
tors, and rj — — 1 (rj — +1) denotes the retarded commu- 
tator (anticommutator) Green's function (see Appendix 
for the detailed definition of Green's function) . The latter 
two (G~~' mn and G^' mn ) are used for relatively com- 
plicated systems, as needed>£ ' 21 i 22 ' 23 i 24 ' 25 In this section, 
according to the previous works, we employ this usual 
set of Green's functions. 

To avoid extra complexity, we consider the case of 
5 = 1/2, where the above usual set is reduced to only 



o -r„ 



-r q /2 r q /2 o 




(10) 



with r q = h + zJm(l — A7 q ) and r q = h + zJm(l — 7 q ), 
where z — 2 (z = 4) is the number of nearest neighbors 
and 7 q = cos q x 

Kj = ^(cosfe + cosq y )} is the Fourier 
factor for a chain (for a square lattice). Here we set the 
lattice constant to be unity. By solving Eq. (jTDJ) , one 
derives the commutator Green's functions 



G 



±mT c 



(11) 



with w q = r q r q . The anticommutator Green's functions 
Gq^ have a pole at u> — and one obtains 



C. 



±- 



lim luG 



±- 



A: 



A 



q, + l 



Then the spectral theorem [Eq. (|A6[) ] gives 



A 



q,+l 



^coth^S 
uj a 2 



(12) 



(13) 
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and 



and 



A q7+1 = - 



(14) 



Note that the applications of the spectral theorem to 
G+~ and G~~ yield the same equation [Eq. (fT3|) ]. More- 
over, the following expression is derived from Eq. (|13D : 



-^2\ 



{(sn 



1 x mP. 



v-!^ coth ^a. (15) 



Here the sum runs over the first Brillouin zone. In addi- 
tion to the relation for general spins S, 

{St, Sr} = 2[S(S+l)-(S?) 2 ], (16) 

the spin-1/2 operators have the following properties: 



(st) 2 = (sn 2 = o, (stf = ]. 



(17) 



If, as usual, we demand that the conditions given by 
Eqs. lfT6|) and (fl7|) are satisfied even within the RPA, 
Eq. (fT5|) becomes 



1 1 x ■> mT. 

2 ~~ N 



-coth^a. 



(18) 



One can now calculate the magnetization m from this 
self-consistent equation. 

The above procedure has been often used in the case 
of systems with an easy-axis anisotropy (A > 1) and 
a transverse magnetic field. 32 Despite this, the obtained 
results are in poor agreement with the numerical results 
as we shall discuss later. 



IV. CONTRADICTION BETWEEN TWO 
EQUATIONS 

In the previous section, according to some pre- 
vious studies, we use the set of Green's functions 
{Gf~ rj ,G~~ ql Gl~ rj \ and adopt the conditions described 
by Eqs. (TlBl) and (jXTJ) . However, this choice of Green's 
functions is not enough to take care of all directions in 
spin space. In this section, we consider all combinations 
of spin operators on the choice of the operators A and B 
in Eq. ((XT]) : 

G il v = (( s ii Sj))r, (M) v = x,y, z). (19) 
The equations of motion are given by 

u>G£» = ([S?,S%] V ) + (([S?,B\;Sj))v,», (20) 

with 



[Sf,H] = -i ^ J ik (S? SI - ASfS z k ), 

k 



(21) 



k 



J ik (ASfS* k -S?S%)+ihS?, (22) 



[Sf , H] = -iY] J ik (SfS x k - S*S%) - ihSf. (23) 



k 



Applying the RPA decoupling expressed by Eq. ([7]) , and 
performing the Fourier transformations given by Eqs. ([8]) 
and ((9|) , one derives the following matrix equation instead 
of Eq. |10]): 



LUl - 







o -*r q o 



sixu 



G yu = k yv 



(24) 



By solving Eq. (|24j) . one obtains the commutator Green's 
functions 



and 



Gzv 
q,-i 



jA v _\ + if q AI«i 



U3* 



(25) 
(26) 

(27) 



Since the commutator of spin operators is local in site 
indices, Aq"_ 1 are independent of q and thus we drop the 
subscript q from A^l-p Here we use the fact A 2 !^ = 0, 
which is valid for any v. Additionally, one obtains 

/~ixv kxv /~<yv stzv n /oa\ 

U q, + 1 — yV q, + l' °q,+l — °q, + l — U 

from the anticommutator Green's functions. Then the 
spectral theorem [Eq. (|A6j) ] gives 



and 



yi q,+i 



A q,+i 



^coth^a. 



-*T q A^ 



coth 



(29) 



(30) 



Incidentally, the application of the spectral theorem to 
Gq 1 ^ gives the trivial equation A q ^ +1 = Aq 1 ^. Then 
Eqs. (|29| and ([30| become, for v = y and v = z, respec- 
tively, 



1 



((sfY) = i- T 



N 2lu c 



(31) 



and 



((S?f) = -J- V^coth 
ww/ 2uj a 



(32) 



If we demand that the property specific to the spin-1/2 
operators, (Sf ) 2 = 1/4, is satisfied even within the RPA, 
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Eq. (fl8|) is derived again from Eq. (|3ip. Additionally, 
from Eq. (|32p with (iSf) 2 = 1/4, one more equation 

1 lyrtf* CQth ^a (33 ) 
q M 

is derived, which cannot be obtained by using only 
{C; r , r Gr.", G?r }. However, the two Eqs. (HHJ) and (02D 
obviously contradict each other, except in the isotropic 
(A = 1) case. In other words, the different values of the 
magnetization m are obtained from each equation. This 
fact means that the above procedure with the setting 
(Sf ) 2 = (Sf) 2 = 1/4, which is known as a good approxi- 
mation for isotropic models, is not valid for a system with 
an easy-plane exchange anisotropy. The same problem 
should arise whenever there is an anisotropy in the plane 
perpendicular to the direction of the magnetization (i.e., 
whenever the system is not rotationally invariant around 
the direction of the magnetization). 

This problem may be avoided by assuming that a cer- 
tain restriction has to be imposed on the choice of Sj 
in Eq. (fT9|) . This idea was introduced by Brown^ for 
the easy-axis case. If the choice of Sj is restricted to 
be perpendicular to the anisotropy axis, i.e., Sj = Sj, 
only Eq. JTHJ is obtained from Eq. (J29|) while Eq. {30]) 
yields the identity = 0. In contrast, under the re- 
striction that = Sj, only Eq. (|3"3"|) is obtained from 
Eq. (|30|) . In Fig. [I] the results for the temperature depen- 
dence of m obtained from each case are compared with 
the QMC results. In the QMC calculations, where we use 
the continuous-time loop algorithm ; 3 ^ 34 ' 35 z direction is 
taken as the quantization axis. The transverse magneti- 
zation in x direction is measured by using the improved 
estimator technique^ For each temperature, measure- 
ment is performed for 6.6 x 10 5 Monte Carlo steps after 
discarding 8 x 10 3 steps for thermalization. The system 
sizes are 128 and 32 x 32 for the chain and the square 
lattice, respectively, which are large enough to produce 
the data in the thermodynamic limit in the resolution of 

Fig. m 

We can see at once that the results obtained from 
Eq. (|3"2")) . which are shown in Fig. [IJb), violates the in- 
equality | (Si) | < S. Thus the assumption that the choice 
of Sj is restricted to be Sj is obviously not correct. On 
the other hand, as seen in Fig.[lja), the choice Sj = Sj 
(or the procedure of the previous section) yields more 
reasonable results in the sense that |(Sj)| < 5 is satis- 
fied. However, in comparison with the QMC results, the 
approximation gets worse as the temperature increases. 
In particular, the intersections of the lines for the differ- 
ent anisotropy parameters are found in the QMC results, 
whereas the lines in Fig. [lja) do not cross each other. 

After all, the results obtained from either choice of Sj 
do not agree with the QMC calculations. In the first 
place, the reason that a restriction on the choice of Sj is 
present is unclear. The result of the comparison between 
Figs. [Ha) and [He) indicates that the procedure of the 
previous section is inappropriate to treat anisotropic spin 




' 0.0 0.5 1.0 1.5 2.0 



temperature T/J 

FIG. 1: The temperature dependence of the magnetization 
of an easy-plane ferromagnet at h/J — 0.1 for A = 1 (the 
isotropic Heisenberg model), A = 0.7, A = 0.4, and A = 
(the XY model). Comparison between the results obtained 
from (a) Eq. (j"8]l. (b) Eq. i[33]l. and (c) the QMC calculations. 
The error bar of the QMC results are smaller than the line 
width. 



systems, although it has been often used for the case 
of systems with an easy- axis anisotropy (A > 1) and a 
transverse magnetic field^ i 21 ' 22 i 23 i 24 i 25 



V. DISCUSSION AND ATTEMPT TO 
IMPROVE THE METHOD 

We now discuss why the contradiction between 
Eqs. (fT8|) and ([33]) occurs, and attempt to improve the 
method. The exact values of ({Sf) 2 ) and ((Sf) 2 ) are 
equivalent (= 1/4) from the properties of the spin-1/2 
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operators. However, since we adopt the decoupling ap- 
proximation, of course, ((Sf) 2 ) and ((Sf ) 2 ) have errors 
from the exact values, respectively. As seen from Eq. (JlJ, 
the parameter A induces the anisotropy in yz plane, 
and the model is not symmetric under the exchange 
S y <-» S z . Thus, within the RPA [Eq. {7}], the errors 
of ((Sf) 2 ) and ((Sf) 2 ) should also be "asymmetric", i.e., 
|((Sf) 2 >R P A - 1/4| ? \((Sf) 2 ) RPA - 1/4=1, which causes 
the contradiction between Eqs. (|31[) and (|32|) with the 
setting ((Sf) 2 ) = ((Sf) 2 ) = 1/4, namely Eqs. iJTSJ and 
(f3"3"| . Here the subscript "RPA" expressly denotes the 
approximate values obtained within the RPA. 
For general spin S, in addition to the constraint 



S 2 = (Sf) 2 + (Sf) 2 + (Sf) 2 = S(S + 1), 

we have the operator identities 

s 

n (^-pho, 

p=-S 

n (5?- P )=o, 

p=-S 

and 

s 

n (^-p)=°- 



(34) 

(35a) 
(35b) 



(35c) 



As is mentioned above, within Green's function theory 
with the RPA, these severe conditions cause the contra- 
diction between Eqs. fT5]) and (13"3")) . Therefore we have 
to think about relaxing the local restrictions on the spin 
operators. As is clear from Eq. (f24|) . the calculation of 
G%j is not necessary to derive Eqs. (f2U)) and ([3"0|) . Thus 
the decoupling approximation to the time evolution of 
Sf [Eq. (|2"Tj) ] is not actually required. Meanwhile, as 
for the time evolutions of Sf and Sf [Eqs. [J22]) and (|23|) ]. 
one needs to employ the decoupling approximation to the 
higher-order terms. This implies that it is inappropriate 
to use Eqs. (135b[) or (|35c[) as a conditional equation. Ac- 
tually, the obtained results are unreasonable as shown in 
Figs.flja) and (b). Thus, in our method, we demand that 
the operator identity in the direction of the magnetiza- 
tion [Eq. ([35a) ] is satisfied. From Eqs. fl34]) and (|35a|) . 
we obtain the following conditional equation: 



((^) 2 )rpa + ((^) 2 )rpa = x 



(36) 



The same condition was used by Aoki in his "self- 
consistent spin-wave approach."— Then, from Eqs. (j3Tj) 
and (|32p with the condition given by Eq. (f3"6")) . the self- 
consistent equation 



J_ \ - m(Tq 

N ^ 
q 



r q ) 



2uj c 



coth^ (37) 



0.5 



0.4 



G 

O 

0.3 

N 

I 0.2 

b o.i 




0.0 
0.5 



0.4 



o 

0.3 

I 0.2 

bo 

a 

B 0.1 



0.0 



(a) 

A = 1 (isotropic) 
A = 0.7 
A = 0.4 
A = (XY) 



chain 




0.0 0.5 1.0 1.5 

temperature T/J 



2.0 



FIG. 2: (a) The temperature dependence of the magnetiza- 
tion obtained from Eq. (|37|l at h/J = 0.1. (b) Comparison 
between the results obtained from Eq. (j3T|> . the QMC calcu- 
lations, and the MFT (at h/J = 0.1 for A = 0.7). 



is obtained. Obviously, when A = 1, this equation re- 
produces the result obtained by the conventional RPA 
theory for the isotropic Heisenberg model as expected. 

As seen in Fig. HJa), the results from Eq. (|37[) are in 
good accord with the QMC results shown in Fig. [TJc). 
Especially, the intersections of the lines for the differ- 
ent anisotropy parameters are found as expected. In 
Fig. H)Jb), we show the comparison between the results 
obtained from Eq. (|3T|) . the QMC calculations, and the 
mean- field theory (MFT). In the MFT, since the terms 
SfSy + ASfS? in Eq. flU are neglected, the obtained 
magnetization curve does not depend on the anisotropy 
parameter: 

m (^ = \t^ h + Z Y MFT)) . 08) 

Meanwhile in our procedure (and in the QMC calcula- 
tions), the effect of the anisotropy can be appropriately 
taken into account. At low temperatures (even at T = 0), 
the magnetization is suppressed by the quantum fluc- 
tuations induced by A. In contrast, at high tempera- 
tures, the parameter A plays an opposite role: since the 
easy-plane anisotropy energetically favors the spin align- 
ment in the easy plane, whereas the quantum fluctua- 
tions have an insignificant effect, the magnetization is 
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enhanced as A decreases. As a result, the intersections 
of the lines for the different anisotropy parameters are 
found in Figs.QJc) and[2^a). 

The application of Green's function theory to Heisen- 
berg ferromagnets with an easy-plane anisotropy was in- 
vestigated also in Refs. [H and[3E In particular, recently 
Hu et alM* calculated the magnetization, susceptibility, 
and transverse correlation functions of the system de- 
scribed by 

H = -\ E J a (S?S*+ASVS]+S?S*)-hJ2 st- (39) 

i,j i 

Obviously, this model is exactly equivalent to Eq. (JT|) un- 
der the rotation of the coordinate system. Using Green's 
functions Gf~„ and G~~, they obtained the same results 

'J s 'I t J-'r ^ 

as those shown in Fig. [Ha) from the condition 

(SiSf) = - - {St). (40) 

Just in the case of the coordinate system described by 
Eq. (|39p . where the magnetization appears along z axis, 
the condition given by Eq. (j4"0"l) produced the same effect 
as Eq. (|36p of our procedure, and they obtained the same 
results even though they use only Gjjl and G~~. How- 
ever, this fact does not mean that this choice of Green's 
functions and condition always yields reasonable results 
for any system. In fact, for the model in the coordi- 
nate system described by Eq. (TJJ), the adoption of the 
severe conditions Eqs. (fTrJ)) and (jTTJ) with the usual set 
{Gf~, G~j~,G*~ v } yields the unreasonable results, as we 
showed. 



VI. SUMMARY 

In this paper, we have investigated the application of 
the many-body Green's function theory with the random- 
phase approximation to the study of easy-plane ferro- 
magnets in an in-plane magnetic field. If there is an 
anisotropy in the plane perpendicular to the direction of 
the magnetization, then special attention is required for 
the choices of Green's functions and the conditions to 
determine the magnetization. First, we calculated the 
temperature dependence of the magnetization for vari- 
ous values of the anisotropy parameter A by using the 
usual set of Green's functions {G^~, G^~, G|T } and 
the conditions {S+,Sf} = 1 and (Sf) 2 = 0. The ob- 
tained results did not agree with the QMC calculations 
on the point that no intersections of the lines for the dif- 
ferent anisotropy parameters were found [see Figs. []Ja) 
and E^c)] , which means that this procedure cannot ap- 
propriately take into account the effect of the anisotropy. 

Next, by considering all combinations of two spin op- 
erators, G^ (/i, v = x,y,z), we derived the additional 
equation, which is not obtained by using the previous 
procedure. Then we showed that the two equations con- 
tradict each other under the conditions (Sf ) 2 = (Sf ) 2 = 



1/4. The same problem should also arise for higher spins, 
as long as the system does not have rotational symme- 
try around the direction of the magnetization. To avoid 
this contradiction within the RPA, we relaxed the restric- 
tions; we demanded that only the operator identity in the 
direction of the magnetization is satisfied. For example, 
in the case that the magnetization appears along the x 
axis, the condition ((SJO^rpa + ((Sf ) 2 )rpa = 1/2 was 
adopted. The results obtained by this procedure were in 
good agreement with the QMC calculations. This means 
that our method can appropriately take into account the 
effect of the anisotropy. 

In this paper, we focused on the spin-1/2 case. Here, 
the extension to higher spins S > 1 will be presented. 
For example, according to our procedure, we have the 
relations Sf = 2 and (Sf + l)Sf{Sf - 1) = for the 
spin-1 case. There is one more unknown, ((Sf) 3 ), than 
S = 1/2, and thus an additional equation is required to 
close the system of equations. To this end, one can adopt 
the following relation: 

({SV, SfSf}) + ({Sf,SfSf}) = 3(Sf) - 2((Sf ) 3 >. (41) 

The expression for the left-hand side is obtained from 
Green's functions {(Sf;S*S»)) and ({Sf;SJSp) after 
some algebra similar to the one given in Sec. HVl In 
the same way, for general spin S, the system of self- 
consistent equations is obtained from the expressions for 
{{Sf,(SfrSf}) + {{Sf,(SfrSf}) (n = 0,1,-.. ,25-1) 
and the conditions given by Eqs. (fM)) and (|35ap . 

Finally, the case of easy-axis anisotropy (A > 
1) is also addressed. In this case, it is well 
known that the transverse magnetic field h in- 
duces the spin reorientation transition for both 
the ferromagnetic (J > 0) and antiferromagnetic 
(J < 0) casesii^^^^^^^^L^^i In or _ 
der to calculate the magnetic properties, many pre- 
vious authors£ i 21 ' 22 i 23 ' 24 i 25 have employed the Green's 
function theory with the set of Green's functions 
{G+~' m ™, G^~' mn , G^ mn } and the severe conditions 
such as Eqs. (fT!)]) and p7[) in the RPA. In other words, 
they have employed the method presented in Sec. IHII or 
its extended versions for higher spins. However, accord- 
ing to the results we showed for the easy-plane case, this 
choice of Green's functions and the conditions is not al- 
ways appropriate. Actually, unacceptable results, espe- 
cially the violation of the property of spin-5 1 systems, 
(Si) | < 5 1 , were found in some previous works (for ex- 
ample, see Ref. [23j; although it is not explicitly stated, 
the same problem is found also in some other works). 
The theoretical treatment of these systems is generally 
more complicated than the easy-plane cases since the 
transverse magnetic field causes changes in the direction 
of the magnetization. Nevertheless, the results reported 
here will be helpful for investigating these problems. 
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APPENDIX A: GREEN'S FUNCTION THEORY 

We briefly review the many-body Green's function the- 
ory. The retarded commutator (rj = — 1) or anticommu- 
tator (rj — +1) Green's function is defined by 



case. Incidentally, Eqs. (|A3|) and (|A4[) can be rewritten 
in the following more convenient form: 



({AB}) 



— lim / dio \G AB (u 



iS) 



(3uo 



-Gif(oj-iS)] coth^+C^f 



(A6) 



and 



([AB]) 



■ lim 



du> 



2ttS^ j _ x 
x [G A f(u; + i5)-G A f(u;-iS)] 



(A7) 



In calculations of spin systems, Eq. (|A6j) generally plays 
the role of the self-consistent equation, whereas Eq. (|A7|) 
is trivial and gives no informations. 



G AB (t-t') = ((A(t-t');B)) n 

= -i6(t-t'){[A(t),B(t%), (Al) 

where 6(t — t 1 ) denotes the step function and (• • • } = 
r Ti{e^^ H ■ ■ ■ )/Tr(e~^ ff ) is an average over the ensemble 
with the inverse temperature (3 — 1/T. [A, £>]-i = [A, B] 
and [A, B] + \ = {A, B} denote the commutator and anti- 
commutator of operators A and B, respectively. Green's 
function satisfies the following equation of motion in en- 
ergy space: 



u>G AB (u) = {[A,B] ri ) + (([A,H};B)) 



(A2) 



To obtain the solution of G ab (oj), a certain decoupling 
approximation of higher-order Green's functions in the 
right-hand side of the above equation is usually adopted. 
Then the expectation values and correlation functions 
can be calculated from the spectral theorem^ 



(AB) = -L, ta /V*f(" + ")-Off(»-") 



AB 



r<AL 



(A3) 



and 



x 1 2tt<5^o 



AB 



c<al 



(A4) 



where 



C AB = lim wG 



AB 



(A5) 



If the anticommutator Green's function has a first-order 
pole at oj = 0, the terms C AB /2 in Eqs. (TJ3| and ([54)) 
are required. Meanwhile, it is well known that the com- 
mutator Green's function has no pole at lu = in any 



8 



* Electronic address: yamamoto@kh.phys.waseda.ac.jp 
^ URL: |http : / /w ww . kh . phys . wase da. ac . jp/ ] 

1 S. V. Tyablikov, Ukr. Mat. Zh. 11, 289 (1959); S. V. 
Tyablikov, Methods in the Quantum Theory of Magnetism 
(Plenum Press, New York, 1967). 

2 H. B. Callen, Phys. Rev. 130, 890 (1963). 

3 R. A. Tahir-Kheli, Phys. Rev. 132, 689 (1963). 

4 J. A. Copeland and H. A. Gersch, Phys. Rev. 143, 236 
(1966). 

5 R. H. Swendsen, Phys. Rev. B 5, 116 (1972). 

6 A. Oguchi, Prog. Theor. Phys. 46, 63 (1971). 

7 P. Frobrich and P. J. Kuntz, Phys. Rep. 432, 223 (2006). 

8 J. Kondo and K. Yamaji, Prog. Theor. Phys. 47, 807 
(1972). 

9 H. Shimahara and S. Takada, J. Phys. Soc. Jpn. 60, 2394 
(1991). 

10 I. Junger, D. Ihle, J. Richter, and A. Kliimper, Phys. Rev. 
B 70, 104419 (2004). 

11 T. N. Antsygina, M. I. Poltavskaya, I. I. Poltavsky, and K. 
A. Chishko, Phys. Rev. B 77, 024407 (2000). 

12 Y. G. Rudoi and Y. A. Tserkovnikov, Teor. Math. Phys. 
14, 75 (1973); 15, 600 (1973); J. M. Wesselinowa and V. 
P. Ivanova, Phys. Rev. B 54, 9846 (1996). 

13 T. L. Reinecke and R. B. Stinchcombe, J. Phys. C: Solid 
State Phys., 7, 527 (1974). 

14 M. E. Fisher and D. R. Nelson, Phys. Rev. Lett. 32, 1350 
(1974). 

15 H. Rohrer, Phys. Rev. Lett. 34, 1638 (1975). 

16 D. P. Landau and K. Binder, Phys. Rev. B 24, 1391 (1981). 

17 K. Ohgushi and Y. Ueda, Phys. Rev. Lett. 95, 217202 
(2005). 

18 M. Holtschneider, W. Selke, and R. Leidl, Phys. Rev. B 
72, 064443 (2005); M. Holtschneider, S. Wessel, and W. 
Selke, ibid. 75, 224417 (2007). 

19 Y. G. Rudoi and Y. A. Tserkovnikov, Teor. Math. Phys. 
25, 1073 (1975). 

20 E. B. Brown, Phys. Rev. B 36, 266 (1987). 

21 P. Frobrich, P. J. Jensen, and P. J. Kuntz, Eur. Phys. J. B 



13, 477 (2000). 

22 P. Frobrich, P. J. Jensen, P. J. Kuntz, and A. Ecker, Eur. 
Phys. J. B 18, 579 (2000). 

23 P. Frobrich and P. J. Kuntz, Eur. Phys. J. B 32, 445 (2003). 

24 P. J. Jensen, S. Knappmann, W. Wulfhekel, and H. P. 
Oepen, Phys. Rev. B 67, 184417 (2003). 

25 H.-Y. Wang, C.-Y. Wang, and E.-G. Wang, Phys. Rev. B 
69, 174431 (2004). 

26 J.-S. Caux, F. H. L. Essler, and U. Low, Phys. Rev. B 68, 
134431 (2003). 

27 D. V. Dmitriev and V. Ya. Krivnov, Phys. Rev. B 70, 
144414 (2004). 

28 S. Schwieger, J. Kienert, and W. Nolting, Phys. Rev. B 

71, 024428 (2005). 

29 M. G. Pini, P. Politi, and R. L. Stamps, Phys. Rev. B 

72, 014454 (2005); S. Tacchi, A. Stollo, M. Madami, G. 
Gubbiotti, G. Carlotti, M. G. Pini, P. Politi, and R. L. 
Stamps, Surface Science 600, 4147 (2006). 

30 P. J. Jensen, K. H. Bennemann, D. K. Morr, and H. 
Dreysse, Phys. Rev. B 73, 144405 (2006). 

31 D. Yamamoto and S. Kurihara, Phys. Rev. B 75, 134520 
(2007). 

32 Concretely speaking, for example, our Eq. (|18[) corre- 
sponds to Eq. (28) of Ref. [23] through the substitution 

h x — > r q , h x — > r q , h z — > 0. 

33 |http://wist aria.co mp-phys.org/alps-looper/ 1 

34 S. Todo and K. Kato, Phys. Rev. Lett. 87, 047203 (2001). 

35 A.F. Albuquerque et al, J. Mag. Mag. Mat. 310, 1187 
(2007). 

36 R. Brower, S. Chandrasekharan, and U. J. Wiese, Physica 
A 261, 520 (1998). 

37 T. Aoki, J. Phys. Soc. Jpn. 61, 2909 (1992). 

38 Y. G. Rudoi and Y. A. Tserkovnikov, Teor. Math. Phys. 
19, 491 (1974). 

39 A.-Y. Hu, Y. Chena, and L.-J. Peng, J. Magn. Magn. 
Mater. 313, 366 (2007). 



